An Accurate Curved Boundary Treatmentin the Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is an alternative kinetic method capable of solving hydrodynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and natural to parallelize. Because the method often uses uniform regular Cartesian lattices in space, curved boundaries are often approximated by a series of stairs that leads to reduction in computational accuracy. In this work, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary. The proposed treatment of the curved boundaries is an improvement of a scheme due to Filippova and Hanel. The proposed treatment for curved boundaries is tested against several ow problems: 2-D channel ows with constant and oscillating pressure gradients for which analytic solutions are known, ow due to an impulsively started wall, lid-driven square cavity ow, and uniform ow over a column of circular cylinders. The second-order accuracy is observed with solid boundary arbitrarily placed between lattice nodes. The proposed boundary condition has well behaved stability characteristics when the relaxation time is close to 1=2, the zero limit of viscosity. The improvement can make a substantial contribution toward simulating practical uid ow problems using the lattice Boltzmann method.